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Abstract 

The calculation of D-wave charmonium prodution rates in B meson deacys 
under the NRQCD factorization formalism is presented. We hnd that inclusion 
of the color-octet contributions permits us to detect the D-wave charmonium 
states in B decays at present experimental facilities. The same amount signals 
of 2 D-wave state as that of xf' could be observed at LEP and CESR. 
We also predict the relative production rates for four D-wave states are 6 ^ : 

: 62 ■ = 2.5 : 3 : 5 : 7, where 5^, 62 , represent respectively the 

2“+, 1 , 2 , 3 states. 
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I. INTRODUCTION 


Studies of heavy quarkonium production in high energy collisions provide important im- 
formation on both perturbative and nonperturbative QCD. In recent years, a rigorous fram- 
work for treating quarkonium production and decays has been advocated by Bodwin, Braaten 
and Lepage in the context of nonrelativistic quantum chromodynamics (NRQCD) [Q. In this 
approach, the production process is factorized into short and long distance parts, while the 
latter is associated with the nonperturbative matrix elements of four-fermion operators. 
This factorization formalism provides a new production mechanism called the color-octet 
mechanism, in which the heavy-quark and antiquark pair is produced at short distance in a 
color-octet conhguration and subsequently evolves nonperturbatively into physical quarko¬ 
nium state. This mechanism is hrst considered to cancel the infrared divergence in the 
calculation of P-wave charmonium production in B decays p|. But the most important 
progress is that by including the color-octet production mechanism one might explain the xjj' 
{J/'ip) surplus measured by CDF at the Tevatron 0. In the past few years, applications of 
the NRQCD factorization formalism to J production at various experimental facilities 
have been studied [Q. 

Recently, a calcultion of P-wave charmonium production in decays suggests a crucial 
test of color-octet production mechanism p. The authors show that the production rates 
due to color-octet gluon fragmentation to P-wave charmonium states are 2 ~ 3 orders larger 
than that due to the dominant color-singlet quark fragmentation |^. They also predict that 
at the Tevatron the same amount of magnitude of direct prodution rate of 2 P-wave 
charmonium state as tjj' could be observed []^. All these progresses show that the color-octet 
mechanism is crucial important to P-wave charmonium production because the color-octet 
contributions are over two orders larger than the color-singlet contributions. While in S- or 
P-wave states production, the situations are not so critical. This huge divergences between 
color-octet and color-singlet contributions in P-wave states production are most helpful in 
distinguishing the two production mechanisms in experiments and will provide a crucial test 
of color-octet production mechanism. On the other hand, even if the color-singlet model 
predicts the P-wave production rates too small to be visible, they could now be detected 
after and only after including the color-octet production mechanism. So it is the color-octet 
mechanism that could make it possible to search for the P-wave heavy quarkonium states. 
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and then will complete the stndies of charmoninm and bottomoninm families. 

In this paper, we will calcnlate D-wave charmoninm prodnction rates in B decays, as 
a complementation to previons stndies ^ |]^, and also as an independent determination of 
the color-octet matrix elements associated with color-octet D-wave charmoninm prodnction. 
Color-singlet S- and P-wave charmoninm prodcution in B decays have been stndied in the 
literatnres |Q]. The color-octet S- and P-wave prodnction have also been investigated [Q 
i m- /^From the following calcnlations, we will show that the color-octet prodnction 
mechanism is also crncial important to P-wave charmoninm prodnction in B decays. The 
rest of the paper is organized as follows. In Set.II, we give the factorization formnlas nnder the 
NRQCD formalism. In Sec.Ill, we estimate the branching fractions of P-wave charmoninm 
prodnction in B decays by assnming that the valnes of color-octet matrix elements are taken 
by the NRQCD velocity scaling rnles. We also disccnss the detection of P-wave charmoninm 
states in B deacys at present experimental facilities in this section. A short snmmary will 
be given in Sec.IV. 


II. FACTORIZATION FORMULA 

According to the NRQCD factorization formalism, the patial width for the inclnsive 
prodnction of a given charmoninm state if in 6 decays has the following factorization form, 

T(b^H + X) = f (6 ^ cc(n) + X) (o«) . (1) 

n 

Here, r„ is the short-distance snbprocess rate for prodncing a cc pair in conhgnration denoted 
by n (inclnding the angnlar momentnm and the color index 1 or 8). is fhe 

long-distance nonpertnrbative matrix element demonstrating the probability of the cc pair 
in n conhgnration evolving into a physical charmoninm state H. While T^ can be calcnlated 
pertnrbatively as an expansion in conpling constant as{Mb), {On^ is a nonpertnrbative 
parameter, and practically can only be determined by htting to the experimental data. 
Fortnnately, nnder the framework of NRQCD, the relative size of {On') can be ronghly 
estimated by nsing its scaling property in (controled by the velocity scaling rnles), where 
V is the typical relative velocity of the heavy qnark inside the bonnd state H. 

The short-distance snbprocesses h —>■ ccg (with q = s,d) are described by the effective 
Hamitonian 
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( 2 ) 


(-1 o/^ _ n 

Heff = -^VcbV *,{—-^C7^(l - 75 )cs7m(1 “ 75)& 

+{C^ + C.)c-f,{l - 75)rW(l - 75)T“6), 


where Gp is the Fermi constant and are KM matrix elements. The coefficiets Cp and 
C- are Wilson coefficients at the scale of /r = To leading order of as{Mb) and to all 
orders of as{Mi,)ln{Mw/Mf,), they are 

C+{Mb) ^ (3) 

C.{Mb) ^ [as{Mb)/as{Mw)Y^^^^. (4) 


In this effective Hamitonian Eq.(^, the hrst term provides the color-singlet contribntion of 
charmoninm prodnction in b decays, he., it prodnces cc pair in a color-singlet state, while the 
second term prodnces cc pair in a color-octet state and provides the color-octet contribntion. 

To D-wave charmoninm prodnction, as argned in |^ |^, the color-singlet and color-cotet 
processes are both scaled with the same orders in In 6 decays, the color-octet contribn- 
tions inclnde snbprocess for spin-singlet and and snbprocesses for spin-triplet 
D-wave states. In the factorization formula Eq.(p, the associated matrix elements are all in 
and they have following relations according to the NRQCD velocity scaling rules. 



(5) 



( 6 ) 


Here, the symbol represents the physical spin-singlet D-wave charmoninm state, and 6 ‘j 
(J = 1, 2, 3) for spin-triplet 1 , 2 , 3 states. The notations represent the cc 

pairs conhgurations with angular momentum L = 2. Because the color-singlet and color- 
octet processes are all in the same orders in at leading order in as and all these 
processes must be taken into account for a consistent calculation. The color-singlet matrix 
elements can be related to the second derivative of the nonrelativistic radial 

wave function at the origin |i?^(0)p for D-wave by 


(OfCD,)) 

(oi^eoj)) 


15(2J+ 1)W 

47r 
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(7) 

( 8 ) 





The wave function may be evaluated by potential model so that the color-singlet matrix 
elements are hxed. Whereas the color-octet matrix elements \ and {O^'’\ are 

free parameters and must be determined by experiment. 

For the spin-singlet D-wave charmonium state, at leading order in ag, the color-singlet 
subprocess can not contribute to its production in b decays. At leading order in the 
production can only come from the color-octet subprocess. Using the effective Hami- 
tonian(ll), one can calculate the short-distance coefficient P], and get 


r(6 ^ ccCS^^^) + s,d^S^ + X) = 


3(orC5o) 


2M2 


(C'+ + 


with 


(9) 


f° = lU 
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Gl 
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4M2 

ur 


( 10 ) 


Here the relation iWsP + iKdP ~ 1 has been used. 

For spin-triplet D-wave states, at leading order in the color-singlet process can only 
contribute to the production of state ^Di. Other two states (^£* 2 , can not be produced 
in color-singlet process and can only be produced in color-octet processes. 

We hrst calculate the color-singlet contribution to ^Di charmonium production in b de¬ 
cays. This color-singlet production rate can be calculated by making use of the covariant 
formalism |p . /^From standard BS wave function in nonrelativistic approximation for vector 
mesons, the bound state wave function can be project out as. 


<?) = '^{JM\lsLm) 4(1 + 

sm 


( 11 ) 


To D-wave function, <F(P, k) must be expanded to the secdond order in the relative momen¬ 
tum q. The first and second derivative of the wave function are 

4a(9) = y {JM|lsLm)[Af„7„ MM+ f>) + M, MM+ f )7dfe™(«); (12) 


ZlVl^lVl 

After integrating over the relative momentum q, the orbit angular momentum part of the 
wave function will depend on the radial wave function or its derivatives at the orgin ^^(O), 
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Rp{0), and -R'^(O), respectively, for S'-wave (L = 0), P-wave (L = 1), and P-wave {L = 2) 
states, 


^ ^ ^oo(q) = ^=^s(0). 


(27r 


■\/4^ 


(14) 




(16) 




(16) 


where the polarization tensor’s label m is magnetic qnantnm nnmber. For the spin-triplet 
case where J = 1, 2, 3, nsing explicit Clebsch-Gordan coefficients, there are following 
relations for three cases [ITT . 




AMR 


AMR 


^\yal3 4]\^2pp h 


sm 

^(3J^|ls2m) 




2 ^M, 


'Ran ^rPpa'P 9 + ^pa ^rapa'P 9 ) ■> 


M)As) ^ AJz) 

^aP ^p ^aPp- 


(17) 

(18) 
(19) 


Here, Cq,, eap, e^pp are spin-one, spin-two and spin-three polarization tensors which obey the 


projection relations [11 


E ,(m) (m.) _ / „ I PaPP^ _ ^ 

^P - \.-9aP + = IRp, 

^^^aP^^^P' = 7:[Paa''Ppp' +VaP'Vpci'] - RPapRa'P', 

m ^ 

^ -p{'Paa'Vpp'Vpp' + 'Paa''Ppp''Ppp' + VaP'RPa'Vpp' 

m 

+ 'PaP'Vpp'Vpa' + 'Pap''Ppp''Ppa' + Vap'RPa'VpP') 

— —{'Pap'Ppa''Pp'p' + 'Pap'PpP'Ra'p' + 'Pap'Ppp'Ra'P' 

+ RapRPa'RP'p' + 'Pap'Rpp'Ra'p' + RapRPp'Ra'P' 

+ 'PppVaa'Vp'p' + 'Ppp'RaP'Ra' p' + 'PppRap'Ra'P') ■ 

Using these relations, the matrix elements < 0\{cc)v-AfDj > can be calcnlated ont. 


( 20 ) 

( 21 ) 


( 22 ) 
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< 0|(cc)y_Ap-Dl > 

< Q\{cc)v-a?D 2 > 

< Q\{cc)v-a\^D 2 , > 


20^3 i?^(0) 

= 0 , 

= 0 . 


(23) 

(24) 

(25) 


Adopting this results, we obtain the partial width for color-singlet 6 ^ production in b decays, 


r{b ^ cc(‘‘£)i‘’) + s.d ^ + A') = 


25 |flj,(0)|% cmA"5(i 

g„2 ^,^5 ">(20+ C_)(l+ )(1 ) (26) 


5(0f(=-Di)) 

9 M6 


(2C+-C_)2(l + ^)fo. 


(27) 


The color-octet contributions to 5} production in b decays are carried out by using the 
same way H] 


r(6 — {cc)s s,d —^ (5)^ -|- X) = ( 

r(6 —> {cc)s s,d —^ ^2 + X) = ( 

r(& —>■ (cc)s + s,d ^ 6^ + X) = ( 






2M2 

- + ■ 

Mt 

(o?(“Si)) 


{ofep,}) 

2M^ 

M} 

{o?es,)) 


(o?ep.)) 


Q 1 \/T 2 

){c^ + c.ni + ^)r^, (28) 

){c^ + c.ni + ^)r^, (29) 
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){c^ + c.ni + ^)r^. (30) 


Adding the color-singlet and color-octet contributions together, for production, the patial 
width is 

r(& + X) = (1 + - c*-)^ 


+( 


of (3^1) 
2M2 


Mf 


){C^ + C.f]. 


(31) 


III. NUMERICAL ESTIMATION AND DISCUSSION OF THE DETECTION 


In the numerical calculations, we take the input parameters as 

Mb^ Mb = 5.3GeV, M = 2Mc = 3.8GeV, a,(Mfe) = 0.20, a,(M^) = 0.116, (32) 

and then. 
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C+{Mb) = 0.87, C-{Mb) = 1.34. (33) 

By using the standard method 0, we obtain the inclusive branching fractions for and 6 j 

?o) 


BR{B ^6^ + X) = 0.76 
BR{B ^S^^+X) = 0M9 
BR{B ^S^ + X) = 1.7( 


Off'S'n 


M3 




ofcs,; 


2M3 


+ 


mF~ 


), 


BR{B ^5l + X) = 2.4( 


2M| 

of(3^i; 




+ 


MJ 

of(’n) 


). 


(34) 

(36) 

(36) 

(37) 


2M| M5 

where the heavy quark spin symmetry relations for the color-octet matrix elements have 
been used, 

2J+ 1 




1^} ('3 I 

Of^(3R 


3 

2J + 1 


ofl^Si)), 


ob’-Pi) 


(38) 

(39) 


/^From the above results, we can see that the color-octet contributions are crucial important to 
D-wave charmonium production in B decays. The color-singlet contribution is not vanishing 
only in the case of J = 1. Furthermore, in that case, the color-singlet contribution is too small 
compared with the color-octet contributions, because the color-octet matrix elements are at 
the same order in as the color-singlet matrix element but the short-distance coefficients 
in the color-octet terms are about 50 times larger than that in the color singlet term. So, 
the color-octet contributions are much more important than the color singlet contribution. 

The values of color-singlet matrix elements are gotten from the potential model calcula¬ 
tion |i?^(0)p = 0.015GeV'^ 1^. Using the relations Eqs.(8-9), we get 

(of (^Ds)) = 0.27GeUf = 0.16GeU^ (40) 

We estimate the size of the color-octet matrix elements by using the naive NRQCD velocity 
scaling rules. 


M| 

M3 




Pi] 


Ml 

Of (3^2) 
M7 


of (3^1) 
MT 


= 0.0018, 


= 0.0030. 


(41) 

(42) 
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Here, these equation only make sense in the estimate of order of magnitudes (we may also 
use another estimate for these relations as in Ref. , in that case, the results of this paper 
will be enhanced about one order). Adopting these matrix elements values, the predicted 
branching fractions for D-wave charmonium states are. 


BR{B - 

^S^ + X)= 0.23%, 

(43) 

BR{B - 

(5)= + X) = 0.28%, 

(44) 

BR{B - 

+ 0.46%, 

(46) 

BR{B - 

+ X) = 0.65%. 

(46) 


We can easily hnd that the relative production rates predicted above are 6 ^ : 51 : 62 '■ = 

2.5 : 3 : 5 : 7. If we do not take into account the color-octet mechanism the relative rates 
would be : (52 : = 0 : 1 : 0 : 0. 

Among the three triplet states of D-wave charmonium, 5^ is the most prominant can¬ 
didate to be discovered hrstly. It is a narrow resonance, and its branching fraction of the 
decay mode is estimated to be p 

B{61 J/V'TT+TT-) ^ 0.12, (47) 


which is only smaller than that of ) by only a factor of 3. Therefore the 

decay mode 5^ —*• J/'i/57r+7r“ could be observable if the production rate of 5^ is of the same 


order as "0'. Comparing the predicted production rate of 2 D-wave charmonium (^) with 
that oi ijj' 0, we can see that the former has the same amount of that of the latter. At 
CESR, the CLEO collobration have observed strong signals of by J/ijjTni triggering 


Recently, at LEP, the OPAL collaboration have also detected rf)' by the same trigger with 
higher statistics [Q. Therefore, we would expect the detection of 2 H-wave charmonium 
state at these two machines will soon be obtained. 

The other two states of the spin-triplet, 61 and 6 ^ are above the open channel threshold 
and are not narrow, and therefore are difficult to detect. However, the analysis of 
spectrum in decays performed by the OPAL collaboration shows some events above the 
background around 3.77GeV besides the events peak at 'ip'{36S6) H- If these events are 
hnally conhrmed to be associated with the 1 H-wave charmonium state -0(3770), they 
may be mostly from b decays. Moreover, 0(3770) could also be seen in a DD (charmed 
meson pair) hnal state. 












At the Tevatron, about 10® b quarks are produced with an accumulated luminosity of 
100P6“^, which implies more than 10® P-wave charmonium particles could be produced. 
The estimated combined branching ratio 5(^2 Ji ~ 7x 10“^ shows 

that more than 10^ events of 2 P-wave charmonium will be detected at this collider. 
Furthermore, the combined branching ratio B{5^ ^Pili ^Pi J/V’tt®, J/il) has 

been estimated to be about lO”"^ [^. This shows that more than 10^ events of this particle 
can also be detected. 


IV. SUMMAY 

We have calculated P-wave charmonium production rates in B decays in this paper. 
We hnd that the inclusion of color-octet production mechanism can help us to obtain the 
P-wave states production rates at an observable level. In our calculations, the values of the 
color-octet matrix elements are taken by assuming that the NRQCD velocity scaling rules 
are valid in these cases. Practically, these matrix elements can be determined by htting 
the theretical prediction rates to the experimental data. In [0], these matrix elements 
also appear in the calculations. If P-wave charmonium states are detected in the future, all 
these matrix elements can be extracted from experimental data, and the compararison of the 
results from different experiments will provide another important test of the universality of 
the color-octet matrix elements, the NRQCD velocity scaling rules, and further the NRQCD 
factorization formalism. 
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